Interstitial bismuth dimers and single atoms 
as possible centres of broadband near-IR 
luminescence in bismuth-doped glasses 



V.O.Sokolov 1 , V.G.Plotnichenko, E.M.Dianov 

Fiber Optics Research Center of the Russian Academy of Sciences 
38 Vavilov Street, Moscow 119333, Russia 



Abstract 

Absorption, luminescence and Raman spectra of interstitial bismuth atoms, Bi°, and 
negatively charged dimers, Bi^ , in alumosilicate, germanosilicate, phosphosilicate 
and phosphogermanate glasses networks are calculated by time-dependent density 
functional method. On grounds of this calculation an extension of our previously 
suggested model of broadband near-IR luminescence in bismuth-doped glasses is 
put forward. 
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Introduction 



Near-IR broadband (1100 - 1400 nm) luminescence in bismuth-doped glasses 
discovered in Refs. [1,2] is being studied intensively. By now the lumines- 
cence has been observed in many bismuth-doped glasses, such as alumosilicate 
(e.g. [1-7]), alumogermanate (e.g. [4,8-10]), alumoborate (e.g. [10-12]), alu- 
mophosphosilicate, alumophosphate, alumophosphoborate [4,12-14], chalco- 
genide [15,16] and in several bismuth-doped crystals (RbPb 2 Cl 5 [17], FAU-type 
zeolites [18], BaB 2 4 [19]). The bismuth-related IR luminescence is used suc- 
cessfully in laser amplification and generation (see, e.g. [20,21]). The present 
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state of research and applications of near-IR luminescence in bismuth-doped 
glasses are examined in the recent reviews, Ref. [22,23]. 

However there is no commonly accepted model of the IR luminescence center. 
Several models are suggested, such as electronic transitions in Bi + [11,13-15], 
Bi 2+ [15] and Bi 5+ [2-4, 8] interstitial ions, in BiO interstitial molecules [6], 
Bi 2 , Bi^" and Bi 2- interstitial dimers [5,9,16,24,25], other bismuth clusters [9], 
Bi0 4 complexes with tetrahedral coordination of the central bismuth ion [26] . 

In our opinion there is good evidence for negative bismuth dimers being the 
luminescence centres (see e.g. [7,16,25]). However, up to now, there are no 
convincing explanation of changes in absorption and luminescence spectra 
with glass composition (see e.g. Refs. [21,23]). In the context of the Bi^ and 
Bi^ - dimers as centres of the IR luminescence it would appear reasonable that 
such changes are caused by Stark shift and splitting together with certain 
rearrangement of dimer states in glass network electric field. The latter does 
obviously depend on glass composition. 

Alternatively, there is an evidence that single Bi atoms and ions do occur as 
well in the bismuth-doped glass host. We studied Raman scattering of Ar laser 
457.9 nm-light in optical fibres with bismuth-doped alumosilicate glass core 
and observed a narrow intensive band in the low-frequency part of the Raman 
spectra, near 110 cm -1 (Ref. [25], Fig. 4, "a" band). Our calculation [25] 
of neutral single bismuth atom, Bi°, and negative single bismuth ion, Bi~, 
incorporated in alumosilicate host in the sixfold ring interstitial sites proved 
both Bi° and Bi~ have three vibrational modes practically not mixed with 
any vibrations of the rings. The frequency of these vibrational modes are 
found to fall in the 75 - 105 and 85 - 130 cm -1 ranges, for Bi and Bi°, 
respectively (see Fig. 5 in Ref. [25]). One mode corresponds to displacement of 
the interstitial bismuth atom along the rings axis and the other two correspond 
to transverse displacement of the atom. All three vibrational modes are Raman 
active, scattering in the transverse vibrations being an order as intensive as 
in the longitudinal ones. Hence the Raman band near 110 cm -1 in bismuth- 
doped alumosilicate glass may be indicative of single bismuth atoms, Bi°, or 
negative ions, Bi~ in sixfold ring interstitial sites of the glass network. Recently 
electronic transitions in the Bi° atom has been supposed to be responsible for 
the near-IR luminescence in bismuth-doped glasses [27]. 

The aim of this work were to verify the assumptions regarding both the net- 
work electric field and the interstitial bismuth atoms by immediate calculation 
of optical properties of interstitial Bi atoms and ions in networks of silica- and 
germania-based glasses. 
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Calculations 



In our previous works [24, 25] we used calculated data for free negatively 
charged Bi 2 dimers to explain absorption and luminescence spectra of bismuth- 
doped glass. Experimental data available for free Bi° atom were used to the 
same purposes in Ref. [27]. However an influence of glass network (interstitial 
electrostatic field, in essential) must be taken into account to discuss optical 
spectra of interstitial dimers and atoms. 

We performed quantum-chemical calculation of interstitial electric field in net- 
work of glasses of several compositions (silica, germania, alumosilicate, ger- 
manosilicate, phosphosilicate, phosphogermanate). Cluster approach was used 
to model the atomic environment of interstitial sites. Two coaxial sixfold rings 
each formed by Si0 4 , Ge0 4 (for germanosilicate glasses with low germania 
content), AIO4 or = P03 tetrahedra or GeC>6 octahedra (for germanosilicate 
glasses with high germania content) were incorporated in proper numbers 
in the clusters. The distance between the ring centres was optimized during 
the calculations. To ensure aluminum atoms being fourfold coordinated, extra 
electrons were placed in the aluminum-containing clusters, one electron per 
aluminum atom. In GeOe-containing clusters two extra electrons were added 
per sixfold coordinated germanium atom to make it be stable in this coordina- 
tion. Dangling bonds of the outer oxygen atoms in the clusters were saturated 
with hydrogen atoms. All calculations were carried out with GAMESS (US) 
quantum-chemical code [28] by DFT method using BLYP functional which is 
known to provide nice agreement of calculated geometrical parameters with 
experimental data. We used the bases and effective core potentials developed 
in Ref. [29]. One d-type polarization function with ( = 0.8000 a.u. was added 
in the basis for each oxygen atom. Standard 3-21G basis was used for hydrogen 
atoms. 

Quantum-chemical modelling of configuration and vibrational properties of 
Bi2, Bi^ and Bi^ - dimers and of single Bi° atom and Bi~ ion in alumosilicate 
network was performed in our previous works [24,25] using the above-described 
cluster approach. Calculated configuration for Bi^ interstitial dimer was pre- 
sented in figure 3 of Ref. [25] for alumosilicate glass. In the present work such 
calculations were repeated for silica, germania, germanosilicate, phosphosili- 
cate, and phosphogermanate clusters. The equilibrium configurations of neg- 
atively charged dimers, both Bi^ and Bi^", in the interstitial site formed by 
two sixfold rings are found to be quite similar in all these networks: the dimers 
are aligned along the common axis of the rings, the first bismuth atom being 
practically in the center of one ring and the second atom being between two 
rings. Again, equilibrium position of Bi° atom and Bi~ ion are found to be 
between two six-rings in all these networks. As a representative example, cal- 
culated configuration of alumosilicate glass cluster with Bi° interstitial atom 
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is shown in figure 1. Such configurations turns out to be highly stable for both 
dimers and atom or ion: bismuth atoms do not form bond with any atom of 
the rings and returned to the equilibrium positions even after dimer or single 
atom is displaced considerably from that position. Vibrational properties of 
interstitial Bi° atom and Bi^ dimer in all these networks are found to be quite 
similar to those calculated in Ref. [25] for alumosilicate host (see above). 




Figure 1. Calculated configuration of Bi atom in the sixfold-rings interstitial site 
of alumosilicate glass 

Electric field strength was calculated in the point of the Bi° atom position (or, 
what is the same, near the bismuth atom of the Bi2 dimer which is situated 
between two sixfold rings) and in several (8 - 10) points in its vicinity. The 
value averaged over all such points was used in further calculations of optical 
spectra. The mean electric field strengths in ring interstitial sites calculated 
for several glass compositions are given in table 1. As evident from these data, 
the interstitial electric field is found to be approximately the same in silica, 
germanate and germanosilicate (with low germania content) hosts but consid- 
erably higher both in glasses containing aluminate or phosphate component 
and in germanosilicate glasses with high germania content. This is attributable 
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to charged (A10 4 and Ge0 6 ) or strongly dipole-polarized (0 = P0 3 ) structural 
units in the rings around interstitial bismuth sites in such glasses. 

Table 1 

Calculated electric field in sixfold ring interstitial sites of glass network 



T^ini? comnosition 

X \...L_Li.w wvylll l_/\_j^_j± Ulvll 


Electric field 001 a u 


6Si0 4 


2 


6Ge0 4 


2 


4Si0 4 -2Ge0 4 


2 


4Si0 4 -2A10 4 


5 


4Si0 4 -2Ge0 6 


7 


4Si0 4 -20 = P0 3 


7 


4Ge0 4 -20 = P0 3 


10 



Optical spectra calculations were performed by time-dependent density func- 
tional theory (TDDFT) method using Octopus program [30] and Hartwigsen- 
Goedecker-Hutter pseudopotentials [31] with spin polarization and spin-orbit 
interaction taken into account. PBE density functional [32] was used in the 
ground-state calculation. To obtain the linear optical absorption spectrum of 
the system, the Octopus code excites all frequencies of the system by giving 
certain (small enough) momentum to the electrons and then evolves the time- 
dependent Kohn-Sham equations in real space for a certain real time [33]. 
The dipole-strength function (or the photo-absorption cross section) is then 
obtained by a Fourier transform of the time-dependent dipole moment. Adia- 
batic LDA approximation is used in these calculations to describe exchange- 
correlation effects. The Octopus code uses real-space uniform grid inside the 
sum of spheres around each atom of the system (a single one in our case). The 
sphere radius and the grid spacing were taken to be was 8.0 and 0.25 A, re- 
spectively, in our calculations. The real-time propagation was performed with 
2 • 10 4 time steps with the total simulation time of about 20 fs. The Fourier 
transform was performed using third-order polynomial damping (see [30] for 
details of the code). 

Shown in figure 2 are the calculated optical spectra of free Bi^~ dimer and of 
interstitial Bi^ dimers in sixfold ring interstitial sites of networks of silica (or 
germania or germanosilicate), alumosilicate, phosphosilicate and phosphoger- 
manate glasses. Figures 3 and 4 display similar spectra for free and interstitial 
Bi° atoms and for free and interstitial Bi~ ion, respectively. Compairing the 
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results for free BiJT dimers with previous configuration interaction calcula- 
tions [24,25,34] and those for free Bi° with experimental and theoretical data 
available [35-38] (see figure 5), one readily sees that both transition wave- 
lengths and transition intensities are described adequately in our TDDFT 
calculations. 
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Figure 2. Optical cross-section spectra of free Bi^ - dimer (dots) and of interstitial 
Bi^~ dimers in silica, germania or germanosilicate (with low germania content) host 
(short dashes), alumosilicate host (dashes), in phosphosilicate or germanosilicate 
(with high germania content) host (dashes plus dots) and in phosphogermanate 
host (solid) 

Figure 2 shows calculated optical spectra of interstitial Bi^~ dimer in compar- 
ison with those of a free dimer. Similar spectra for interstitial Bi° atom are 
presented in figure 3. From figure 2 it will be noticed how optical spectrum 
of an interstitial Bi^ dimer is modified in comparison with free dimer. Ac- 
cording to our previous calculations performed by configuration interaction 
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(CI) method [25], absorption at 860, 720, 460 nm and < 400 nm wavelengths 
in the free Bi^ dimer is caused by transitions from the ground state to the 
excited ones. The IR luminescence at 1450, 1300 and 1050 nm wavelengths 
corresponds to spin-forbidden transitions from three lowest excited states to 
the ground state and the visible luminescence at 750 nm is caused to a spin- 
allowed transition with low oscillator strength from one of the other excited 
states to the ground state. Notice that compairing results of CI calculations 
with optical spectra calculated in TDDFT with Octopus package (such as 
shown in figures 2, 3, 4) one should realize that the latters include all possible 
transitions in the system under consideration. 
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Energy, eV 

Figure 3. Optical cross-section spectra of free Bi° atom (dots; the part for energy 
< 1.6 eV multiplied by 10 for clarity) and of interstitial Bi° atom in silica, germania 
or germanosilicate (with low germania content) host (short dashes), alumosilicate 
host (dashes), in phosphosilicate or germanosilicate (with high germania content) 
host (dashes plus dots) and in phosphogermanate host (solid) 
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Figure 4. Optical cross-section spectra of free Bi~ ion (dotted line) and of interstitial 
Bi~ ion in alumosilicate host (solid line). 

For the interstitial Bi^~ dimers, the absorption in the < 400 nm range and 
transitions corresponding to the absorption and visible luminescence in the 
600 - 750 nm is found in our calculations to be practically unaltered in all the 
surrounding glass network compositions studied. The absorption band near 
500 nm is found to be shifted towards shorter wavelengths (near 440 nm in 
phosphogermanate environment) and somewhat intensified. Besides, absorp- 
tion near 800 nm growths. In phosphate-bearing and in germanosilicate (with 
high GeC>2 content) glasses absorption in the 900 - 1000 nm range becomes 
significant. In glasses with no phosphate component present and in germanosil- 
icate glasses with low GeC>2 content, the IR luminescence is not too different 
from that in free dimers and occurring mainly in the 1000 - 1300 nm range. 
When passing to phosphate-bearing compositions and for high GeC>2 content 
one finds the IR luminescence bands to be shifted towards longer wavelengths, 
in the 1400 - 1600 nm range. 
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Figure 5. States, transition types and lifetimes in Bi° atom corresponding to the 
6s 2 6p 3 configuration [35-38]: (a) 876 nm, Ml, 32 ms; (b) 648 nm, Ml, 156 ms; 
(c) 462 nm, Ml, 18 ms and E2, 161 ms; (d) 302 nm, Ml, 137 ms; (e) 976 nm, Ml, 
833 ms and E2, 222 ms; (f) 460 nm, Ml, 8 ms; (g) 564 nm, Ml, 43 ms 

Again, optical spectra of Bi° interstitial atoms are found in our calculations to 
change compared with those of free atoms (figure 3). This change is relatively 
low in alumosilicate glasses where Bi° interstitial atoms give rise to absorption 
bands at < 400 nm, near 500 nm and in 600 - 700 and 900 - 1100 nm ranges. 
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This absorption causes IR luminescence in 1000 - 1300 and 1400 - 1600 nm 
wavelength ranges. However the change turn out to be much more pronounced 
in phosphosilicate glasses and in germanosilicate glass with high germania 
content where the absorption should occur in 500 - 600 nm wavelength range 
and, especially, in 700 - 1000 nm range. The corresponding IR luminescence 
in these glasses is found to arise in 1000 - 1400 nm range and near 2000 nm. 

The above-described modification of optical spectra of bismuth interstitial 
dimers and atoms with composition of surrounding glass network gives an 
insight into origin of variations of absorption and IR luminescence bands dis- 
covered in bismuth-doped glasses [21,23]. In particular, characteristic IR lumi- 
nescence band shift from ~ 1100 nm into the 1200 - 1400 nm range is observed 
while turning from alumosilicate compositions to phosphate-bearing ones. 

Ionization potential of the negative bismuth ion, BP, is known to be rather 
low: e.g. in Ref. [38] the electron affinity of 209 Bi is found to be 7600 cm -1 
or 0.94 eV. Because of this, all transitions in < 1300 nm wavelength range 
should, in general, result in BP ionization with neutral Bi° atom formed and 
free electron occurred in conduction band. Since optical cross-sections in BP 
ion are at least an order as high as those in Bi° atom, one would expect the 
interstitial BP ions are destroyed rapidly under both visible and IR irradia- 
tion. Thus interstitial BP ions in glass may promote formation of interstitial 
negatively charged Bi 2 dimers [25]. On the other hand, recombination lumi- 
nescence in 1500 - 1700 nm range is not improbable owing to conduction 
electrons capturing in Bi° interstitial atoms with BP ions formed. 



Discussion 

We have established, on modelling grounds, that interstitial negatively charged 
bismuth dimers, Bi^ and Bi 2 , - , can actually occur in glass network and that 
absorption and luminescence spectra of the dimers correlate well with the 
experimental spectra of bismuth-doped glasses. On the other hand, bismuth 
interstitial atoms, Bi°, and negative ions, BP, are demonstrated to occur in 
glass network as well. In particular, certain features observed in Raman spectra 
of bismuth-doped alumosilicate glasses are attributable to both Bi 2 dimers and 
Bi atoms (ions) present in glass together. 

So it would appear reasonable that both Bi° interstitial atoms and Bi 2 (Bi 2- ) 
negatively charged dimers occur in in bismuth-doped glasses as two types of 
the IR luminescence centres. Relation between interstitial Bi 2 dimers and Bi 
atoms (and ions) concentrations in the given glass is likely to depend predomi- 
nantly on bismuth dopant content. It should be stressed that in bismuth-doped 
glasses considerable part of the dopant bismuth atoms are likely to be bonded 
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with oxygen atoms of surrounding network. However such bismuth centres 
bear no relation to IR luminescence. 

It follows from the relation between calculated intensities of transitions corre- 
sponding to absorption and IR luminescence in Bi° atom that the IR lumines- 
cence time constant is considerably less than in Bi^ or Bi^" dimers and that it 
decreases rapidly with electrostatic field growth (figure 3). Similar calculations 
of Bi2 dimers in electrostatic field prove the transition intensities to change 
only slightly and the transition wavelengths to change almost as rapidly as in 
Bi° atom with the field strength. Therefore, interstitial Bi° atoms should give 
rise to much faster IR luminescence in comparison with Bi^~ and Bi^ dimers. 
Unfortunately the transition intensities in atoms and dimers are hardly to be 
compared directly in our calculational approach. Rough estimate results in 
ratio of the order of 10. Thus, interstitial Bi° atoms may be responsible for 
"fast" component observed in IR luminescence. 

The calculations described in the present article and in Ref. [24, 25] are con- 
cerned with sixfold ring interstitial sites. The sixfold ring are known to be 
the most abundant ring structure in silica- and germania-based glasses net- 
work. However there are other rings as well in such glasses, mainly five-, seven- 
and eightfold ones. We have performed some calculations for interstitial sites 
formed by pairs of such rings in alumosilicate network to examine equilibrium 
states of bismuth atoms and dimers and electric field strength. It is felt that 
the results and conclusions obtained for the sixfold rings still remain valid for 
sevenfold and, supposedly, in eightfold rings. The fivefold rings turn out to be 
too narrow to hold the interstitial dimers and, supposedly, atoms. 



Conclusion 

In summary, we calculated optical spectra of interstitial negatively charged 
bismuth dimers, Bi^~, neutral atoms, Bi°, and negative ions, Bi~, in network of 
silica- and germania-based glasses with various composition. Considering our 
previous works [24,25], the results of the present calculations allow definite 
conclusions, as follows. 

- Both interstitial Bi^ and Bi^ - dimers and interstitial Bi° atoms and Bi 
ions can occur in bismuth-doped silica- and germania-based glasses. IR 
luminescence in 1000 - 1600 nm in these glasses can be caused by all such 
interstitial centres. 

- IR luminescence decay time in Bi° atoms is estimated to be an order as 
short as in Bi^ or Bi^ - dimers . Hence Bi^~ (Bi^ - ) dimers and Bi° atoms 
can be assumed to cause "slow" and "fast" parts of the IR luminescence, 
respectively. 
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- Interstitial electric field strength dependence of absorption and IR lumi- 
nescence in the interstitial Bi^ (Bi^ - ) dimers is found to differ consid- 
erably from that in the interstitial Bi° atoms. Pronounced dependence 
of the IR luminescence spectra and its excitation spectra on glass com- 
position and bismuth content observed in several experiments may be 
attributable to the interstitial field variations. 

- interstitial Bi° atom is found to capture under certain conditions an elec- 
tron from the glass conduction band with Bi~ formed. Recombination 
luminescence in the 1500 - 1700 nm is likely to occur owing to such cap- 
turing. On the other hand, the interstitial negative Bi~ ions are readily 
destroyed under < 1300 nm irradiation. 

Thus our results of the present work and of Ref. [24, 25] concerning intersti- 
tial bismuth atoms and negatively charged bismuth dimers in bismuth-doped 
glasses provide support for the conclusion that both Bi° and Bi^ (Bi^") occur 
in these glasses as two types of the IR luminescence centres. 
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